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We consider here the massive Thirring model regularized 
with the XY Z spin chain. We numerically calculate the mass 
ratios of particles which lie in the discrete part of the spec- 
trum and obtain results in accordance with the DHN formula 
and in disagreement with recent calculations in the litera- 
ture based on the numerical Bethe ansatz and infinite mo- 
mentum frame methods. We also analyze the short distance 
behavior of these states and evaluate the conformal dimen- 
sions. This paper, taken together with the previous one for 
the sine-Gordon model, confirms the duality relation between 
two models formulated by Klassen and Melzer [Int. J. Mod. 
Phys. A 8, 4131 (1993)]. 

PACS number(s): ll.lO.Kk, 11.25.Hf, 11.15.Tk 



I. INTRODUCTION 

The massive Thirring model (MTM) and sine-Gordon 
model (SGM) are important as a testing laboratory for 
understanding ideas proposed for other more complicated 
field theories. 

In this paper we propose to calculate certain physical 
quantities for the MTM by performing an explicit diago- 
nalization of its lattice regularization with the XY Z spin 
chain. As the first task we want to calculate the masses 
of breathers. The previous calculations have been based 
on the semiclassical method 0], on factorized scattering 
theory or on the Bethe ansatz method [^-^. 

The additional interest in avoiding the previously men- 
tioned assumptions is due to recent criticism ||lC)|-[l2|] [the 
authors claim that there is only one breather in the 
whole attractive region, and with different mass than the 
Dashen-Hasslacher-Neveu (DHN) formula predicts] . The 
same authors challange also the well-known duality rela- 
tion between the MTM and SGM (l|-|l|]. The precise 
meaning and extent of this equivalence was formulated 
by Klassen and Melzer (notice that the models are 
not equivalent when they have a finite size in space). 

One important criticism relates to the use of the so- 
called string conjecture. Indeed, violations of this conjec- 
ture are observed in the literature |jl^. Despite the fact 
that, at least untill now, it was not known that these vi- 
olations affect any relevant results, it would be desirable 
to have a calculation which does not rely on the string 
conjecture. 

It is for this reason that we want to treat the MTM 
without using the above-mentioned assumptions. Our 
approach will be based on direct numerical diagonaliza- 
tion of the XY Z spin chain which is a lattice regulariza- 



tion of the MTM |^,^ . This method is suited for analyses 
of low discrete states in the spectrum, but becomes less 
and less effective when we go to higher states. Such an 
approach was used in the literature for other problems, 
e.g., conformal unitary models perturbed by some rele- 
vant operator p^-|2l[]. 

We also intend to calculate conformal dimensions of 
operators creating breather states. There are conjectured 
values for them pTj ]. By explicit calculation we confirm 
this conjecture for the first breather but get different re- 
sults for the second breather. 

Recently we have performed a similar calculation 
for the SGM. The regularization in this case was the 
XXZ spin chain in a transverse field. The results on 
the masses of breathers and conformal dimensions agree 
as statements on relation of two models would suggest, 
so it maybe considered also as an independent check of 
the SGM-MTM correspondence 



II. MTM AS A MASSIVE PERTURBATION OF 
THE GAUSSIAN MODEL 



The MTM is a (1 + l)-dimensional field theory of a 
Dirac spinor field "0, defined classically by the Lagrangian 

Cmtm = i^T^a^V- - AV^V- - U^if,^l^){4>i^^). (2.1) 



Here A is a dimensionful parameter which sets the mass 
scale in a theory which is conformaly inv ariant when 
A = 0. However, although A enters Eq. ( ^.l| ) as a (bare) 
mass, its mass dimension d\ is not equal to 1 but is deter- 
mined from the (nontrivial) anomalous dimension of the 
field ipij). The dimensionless coupling constant g is scale 
invariant (vanishing beta function), but it is not uniquely 
defined due to the existence of different regularizations 
of the (conserved) current ip^'^tp. Correspondingly, there 
is at least a one- parameter family of definitions of g. Our 
definition will be the same as the one used by Coleman 
|l3| (Schwinger definition). We shall find it more conve- 
nient to use the parameter /3 related to g with 
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(2.2) 



Here (3 is the dimensionless coupling constant from the 
dually related SGM. 

In [|l^ it was shown that the MTM can be viewed as a 
perturbed co nfor mal field theory (CFT) when the second 
term in Eq. (2d) is treated as a (massive) perturbation. 
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We will now repeat here some results of their analyses 
relevant for our discussion. 

An unperturbed theory A = (approached in the UV 
limit) is the Thirring model which is a CFT with central 
charge c 1 and an operator algebra generated by 

Lf = {Vjn,n\m e 2Z,n e Z or m e 2Z + l,n e Z + 1/2} 

(2.3) 

where Vm,nix) are primary fields with conformal dimen- 
sions 



N 



'27r /m/32 \^ 27r / m^^ " 



/32 \^ Stt 



/32 V Stt 



(2.4) 



where we used duality relation (^j). From Eq. ( p.4[ ) we 
can read off the scaling dimensions and (Lorentz) spin of 



HxYZ = Hxxz - CT^^+i) , (3.1) 



n=l 



where cr^'W'^ are Pauli matrices, TV is an even integer, and 
Hxxz is the Hamiltonian of the XXZ spin chain: 



N 



Hxxz = E + «+i + ' (3-2) 



where — 1 < A < 1 [we also use standard p ara metriz ation 
A = -cos 7, so 7 e (0,7r)]. In Eqs. (U) and (jsj) 
sector-dependent B.C.'s should be used: 
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N 



(3.3) 



(3.4) 



■^m.n — '—^m.n \ '—^m.n 



167r /32 



(2.5) 



A whole operator algebra is generated by a quartet of 
fields V±i^±i/2, which are connected to the fundamental 
spinor fields '0, '0 by 
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(2.6) 



Now one supposes that Hilbert space of the full (per- 
turbed) theory is isomorphic to that of the unperturbed 
one. From operator product algebra (OPA) it follows 
that the (properly normalized) perturbing operator in 
the MTM (O) is 



(2.7) 



^V = V^i,^^ = ^(^2,0+^-2,0) 



which means that A has mass dimension d\ — 2 — o?2.o = 
2 — 0^/^71. From the condition of relevancy of the pertur- 
bation, i.e., d\ > 0, we obtain Co lema n's bound 0^ < Stt 
(g > -7r/2). Also, from Eqs. and ( p7| ) wc can 

see that the MTM has a U{1) x Z2 x Z2 internal sym- 
metry group. The U{1) acts as Vm.n ^ 
Z2 and Z2 are generated by R : Vm,n 
R ■ Vm,n Vm-n, respectively. 
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From results of Rcf. it |g3|] follows that the XXZ chain 
( |3.2D with B.C. (3.3) gives in the continuum limit a CFT 
with a space of states equal to that of L/ in Eq. (p^), 
where 



(3 = v/8(7r - 7) . 



(3.5) 



That leads us to the conjecture that, aside from irrelevant 
corrections, 
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in the continuum limit. 

Now, the continuum limit is obtained letting N ^ 00 
and /i — > 0, but at the same time keeping fixed the scaling 
parameter fl: 



(3.7) 



In this limit, the mass gaps of the XY Z chain are ex- 
pected to satisfy a scaling law 



m, = /ii/''^G,(7,/i) = h-'^^G,{-f,p) ■ 



(3.8) 



The scaling parameter jl is connected to L (space exten- 
sion of continuum theory, i.e., MTM). For our purposes 
it is enough to know that /i —> cx) (//—;■ 0) corresponds 
to L — *■ 00 (i 0), respectively. 



III. SPIN CHAIN REGULARIZATION OF THE 
MTM 

It was argued a while ago [^j^ that the MTM on a 
cylinder with proper (antiperiodic) boundary conditions 
(B.C.'s) possesses spin chain regularization given by the 
XY Z spin chain defined by the Hamiltonian 



IV. MASS SPECTRUM 

Our goal here is to calculate the mass ratios of particles 
in the MTM in the L li mit using the connection 
with the XYZ spin chain (3.1). First we must numeri- 
cally calculate the mass gaps of the spin chain for finite 
N and h. Then we must make a continuum limit, i.e., 
take N 00 and h ~f 0, keeping p, fixed. Finally we 
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should make a L — > oo, i.e., jl oo, limit. In prac- 
tice, it is preferable to do the following |19|-|2l[]: first take 
N —>■ oo with h fixed and afterwards extrapolate h —> 0. 
A difference is that in the latter case one does /i — > oo 
before h ^ 0. These limits are performed using the BST 
extrapolation method |24 2£ 



We numerically diagonalized Hamiltonian (3.1) for up 
to 16 sites using the Lanczos algorithm. We are inter- 
ested in the masses, so we only need the zero-momentum 
sector. We should note here that in Ref. |^ it was shown 
that true space translations are generated not by an or- 
dinary translation operator on the spin chain, but by its 
square. From this and the fact that Hamiltonian (3.1) 
commutes with the operator C, Eq. (3.4), it follows that 
we can break the Hamiltonian in the momentum-zero sec- 
tor into four sectors named 0^, tt^, where 0, tt is macro- 
scopic momentum and ± denotes eigenvalue of (— 1)^/^C 
(which can only be ±1 because = 1). We considered 
a number of values of coupling constant in the attrac- 
tive regime {g > 0, i.e., A > 0). The structure of the 
spectrum is in agreement with the DHN prediction, i.e., 
we obtain: vacuum, first breather (Bl), second breather 
{B2) (when it exists), and "continuum" in 0^; fermion 
(F) and "continuum" in 0~; antifermion (F) and "con- 
tinuum" in vr^; "continuum" starting with FF and FF 
in 7r+. Names for the particle states and FF, FF contin- 
uum will be confirmed by results for the mass ratios. But 
even we couldn't make an extrapolation (because of the 
poor scaling in the jl —^ oo limit) for the lowest "contin- 
uum" state in 0+ for values of g where the DHN formula 
predicts that it should be of Bl Bl type, its scaling law 
in the fl ^ limit clearly shows that its scaling dimen- 
sion is the one we expect for the Bl Bl lowest continuum 
state, i.e., di^. We should mention also that spectra in 
0~ and TT~ are exactly degenerated which means that the 
F and F mass gaps are equal even on the lattice, which 
was not the case in similar analyses of the SGM in ||2^ . 

In Figs we present numerical results for the scaled 
gaps Gi for four states: fermion (F), first breather {Bl), 
second breather {B2), and lowest state in the FF contin- 
uum (C) . This is of course a check of the scaling relation 
(3.8). Finally, partially extrapolated mass ratios 



fa (A, h) — lim — — — lim , 

— ,"r°°. Gf 



h fixed 



h fixed 



a e {B1,B2,C} 
(4.1) 



and fully extrapolated mass ratios of the first breather 
fBi(A) - limf,(A,/i), (4.2) 

are given in Tables. together with the DHN predic- 

tions 101 and predictions of Fujita et al. |p^ , pT| . Final ex- 
trapolation ft, ^ was possible only for the first breather 
because scaling of the second breather is worse and asks 
for a larger N (probably N > 24). One can see that our 
results strongly confirm DHN and reject Fujita et al. 



V. UV (CONFORMAL) LIMIT OF PARTICLE 
STATES 

Let us now turn our attention to the opposite, i.e., UV, 
limit of our results for the XYZ spin chain. We men- 
tioned in Sec. Ill that it obtained when /i — > 0. From 
conformal perturbation theory we expect the scaling re- 
lation 



(5.1) 



where da is the scaling dimension of the state a, and Q is 
a well-known normalization factor. 
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From Eq. ( p. If) we can obtain the scaling dimensions 
of the particle states F, Bl, and B2 from the condition 
that Ha should be less singular than Ga- Our results are 
given in Table 0. They differ from those conjectured in 
|17t only for the second breather, which has scaling di- 
mension equal to that of the first breather. These results 
are in agreement with those in |Q for the SGM. This 
result for the scaling dimension of the second breather 
was analyticaly confirmed in [ p6[ using an extension of 
nonlinear integral equation method (NLIE). In Figs. ^-|| 
we show the numeric results for reduced scaling func- 
tions, where we used values from Table ^ for the scaling 
dimensions. 



VI. CONCLUSION 

We have calculated in this paper the masses of breather 
states and the anomalous dimensions of related operators 
for the MTM using spin chain regularization. This is a 
direct numerical calculation independent of assumptions 
such as the semiclassical approximation [|l|, factorized 
scattering theory ||^ , or Bethe ansatz method [H^,^ . On 
the other hand, in a series of papers based on numeri- 
cal calculation within the Bethe ansatz metho d |llf| or 
using the infinite momentum frame technique |1C|], dif- 
ferent results have been claimed. Our calculation con- 
firms the conventional spectrum. In addition we cal- 
culate the anomalous dimensions of operators creating 
breather states. It agrees with conjecture in [0 for the 
first breather but disagrees for the second breather. This 
result is consistent with the previous calculation for the 
sine-Gordon model, i.e., consistent with equivalence re- 
lation between the two models IIitI. 



[1] R. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D 
11, 3424 (1975). 



3 



[2] A. Zamolodchikov and Al. Zamolodchikov, Ann. Phys. 

(N.Y.) 120, 253 (1979). 
[3] H. BergknofT and H. B. Thacker, Phys. Rev. Lett. 42, 

135 (1979). 

[4] V. E. Korepin, Theor. Math. Phys. 41, 953 (1979). 
[5] A. Luther, Phys. Rev. B 14, 2153 (1976). 
[6] M. Liischer, Nucl. Phys. B117, 475 (1976). 
[7] M. F. Weiss and K. D. Schotte, Nucl. Phys. B225, 247 
(1983). 

[8] G. I. Japaridze, A. A. Nersesian, and P. B. Wiegmann, 

Nucl. Phys. B230, 511 (1984). 
[9] C. Destri and H. J. de Vega, Nucl. Phys. B504, 621 

(1997) ; D. Fioravanti, A. Mariottini, E. Quattrini, and F. 
Ravanini, Phys. Lett. B 390, 243 (1997); G. Feverati, F. 
Ravanini, and G. Takacs, Nucl. Phys. B540, 543 (1999). 

[10] T. Fujita and A. Ogura, Prog. Theor. Phys. 89, 23 
(1993). 

[11] T. Fujita, Y. Sekiguchi, and K. Yamamoto, Ann. Phys. 

(N.Y.) 255, 204 (1997). 
[12] T. Fujita and M. Hiramoto, Phys. Rev. D 58, 125019 

(1998) . 

[13] S. Coleman, Phys. Rev. D 11, 2088 (1975). 

[14] S. Mandelstam, Phys. Rev. D 11, 3026 (1975). 

[15] D. Delepine, R. Gonzalez Felipe, and J. Weyers, Phys. 

Lett. B 419, 296 (1998). 
[16] A. Gomez Nicola and D. A. Steer, Nucl. Phys. B549, 409 

(1999) . 

[17] T. R. Klassen and E. Melzer, Int. J. Mod. Phys. A 8, 
4131 (1993). 

[18] F. Woynarovich, J. Phys. A 15, 2985 (1982); O. Babelon, 
H. J. de Vega, and C. M. Viallet, Nucl. Phys. B220, 13 
(1983); C. Destri and J. H. Lowenstein, ibid. B205, 369 
(1982); F. H. L. Essler, V. E. Korepin, and K. Schoutens, 
J. Phys. A 25, 4115 (1992); A. Ilakovac, M. Kolanovic, 
S. Pallua, and P. Prester, Phys. Rev. B 60, 7271 (1999); 
Fiz. B8, 453 (1999). 

[19] M. Henkel and H. Saleur, J. Phys. A 22, L513 (1989). 

[20] M. Henkel, J. Phys. A 24, L133 (1991). 

[21] M. Henkel, J. Phys. A 23, 4369 (1990). 

[22] S. Pallua and P. Prester, Phys. Rev. D 59, 125006 (1999). 

[23] F. C. Alcaraz, M. Baake, U. Grimm, and V. Rittenberg, 
J. Phys. A 21, L117 (1988). 

[24] R. Bulirsch and J. Stoer, Numer. Math. 6, 413 (1964). 

[25] M. Henkel and G. Schiitz, J. Phys. A 21, 2617 (1988). 

[26] G. Feverati, F. Ravanini, and G. Takacs, Nucl. Phys. 
B570, 615 (2000). 



FIG. 3. The same as Fig. | but now for A = 0.75 (or 
= 5.78, g = 3.69). 



FIG. 4. Reduced scaling functions Ha{P, /i) at A = 0.3 (or 
(3'^ — 10.13). The legend is the same as in Fig. |l| 

FIG. 5. The same as Fig. || but now for A = 0.6 (or 
= 7.42). 



FIG. 6. The same as Fig. | but now for A = 0.75 (or 
(5^ = 5.78). 



FIG. 1. Scaling functions Ga[l3,n) for the isolated gaps 
[S, BV) plus two lowest "continuum" gaps (CI, C2) of the 
Hamihonian (|^) at A = 0.3 (or (3^ = 10.13, g = 0.76). For 
this value of the coupling constant the DHN formula predicts 
the existence of one breather. The legend in the upper left 
figure applies to all figures in this article. 



FIG. 2. Scaling functions Ga(/3,pi) for the isolated gaps (5, 
Bl, B2) plus lowest "continuum" gap (C) of the Hamiltonian 
(|j) at A = 0.6 (or (3^ = 7.42, g = 2.18). The DHN formula 
predicts now the existence of two breathers. 
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TABLE I. Estimates for the mass gap ratios fa as a function of /i at A = 0.3 (/3^ = 10.13, g = 0.76). In this regularization 
soUton and antisoUton gaps are exactly degenerated. We also added the DHN prediction (only one breather for this value of 
the coupling constant) and the prediction of Fujita et al. (only one breather for all g > Q). The numbers in parentheses give 
the estimated uncertainty in the last given digit. 



fa 








h 








h^Q DHN 


Fujita 


0.8 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


Bl 


1.6341 (3) 


1.7007(4) 


1.718(1) 


1.730 (3) 


1.734 (8) 


1.74(2) 


1.67 (6) 


1.747(6) 1.745 


1.777 


CI 


1.786 (7) 


2.0013(5) 


2.000 (2) 


2.001 (3) 


1.98(1) 


2.00 (2) 


2.07(5) 


2.000 


2.000 


C2 


1.797(2) 


2.0011 (8) 


1.999 (2) 


2.001 (6) 


2.00(1) 


2.00 (3) 


1.93 (8) 


2.000 


2.000 







TABLE II. The 


same as Table | but 


now for A — 


0.6 {0^ = 


7.42, 5 = 2.18). 






fa 






h 








DHN 


Fujita 




0.8 


0.6 0.5 


0.4 


0.3 


0.2 


0.1 '^^o 


Bl 
B2 
C 


1.187 (6) 

1.2(1) 

1.29(1) 


1.2443 (6) 1.2587 (3) 

1.694(2) 
1.536(4) 1.734(2) 


1.26491105 (3) 
1.807020 (8) 
1.99998 (2) 


1.2638 (5) 
1.8753 (8) 
2.003 (2) 


1.254 (2) 
1.913(4) 
2.00(1) 


1.240(9) 1.24(2) 
1.89 (2) 
1.99 (4) 


1.223 
1.935 
2.000 


1.337 
2.000 
2.000 



TABLE HI. The same as Table | but now for A = 0.75 {P^ = 5.78, g = 3.69). 



fa 


0.8 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


DHN 


Bl 




1.027 (4) 


1.030(1) 


1.0255 (4) 


1.0112 (2) 


0.9870 (4) 


0.948(2) 0.91 (2) 


0.905 


B2 


1.0(4) 


1.21 (8) 


1.35 (2) 


1.485 (6) 


1.5716 (3) 


1.6219 (9) 


1.641 (7) 


1.614 


C 


1.528 (3) 


1.25(1) 


1.360 (4) 


1.553 (4) 


1.803 (2) 


2.005 (4) 


1.97(2) 


2.000 



TABLE IV. The same as Table | but now for A = 0.9 (/3^ = 3.61, g = 7.79). 



fa 










h 












DHN 


Fujita 




0.6 


0.5 


0.4 


0.35 


0.3 


0.25 


0.2 


0.15 


0.1 


Bl 
B2 

C 


0.82(1) 

0.9(2) 

0.99(1) 


0.821 (8) 
1.0(2) 
1.010 (8) 


0.795 (2) 
1.13 (6) 
1.166 (8) 


0.779 (2) 

1.18(4) 

1.218(8) 


0.758(1) 
1.18(1) 
1.285 (8) 


0.7356 (5) 
1.202 (8) 
1.366 (6) 


0.7083 (4) 
1.1964 (9) 
1.487(8) 


0.6763 (4) 
1.187 (2) 
1.70(1) 


0.63250 (7) 
1.163 (2) 
1.987(8) 


0.52 (4) 


0.521 
1.005 
2.000 


0.668 
1.336 
2.000 



TABLE V. Scaling dimensions of particle states in the 
MTM as conjectured from our numerical results. 



State 


Operator 


Scahng 


dimension 


Fermion 


Vl,l/2 


16tv 


■ + 




Antifermion 


Vl,_l/2 




■ + 


TT 


First breather 


^2,0 




il 




Second breather 


^2,0 




£. 
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'i/'ii/'2 continuum 
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o h=0.5 
A h=0.4 
□ h=0.3 
T h=0.2 
o h=0.1 

• h=0.08 
A h=0.03 
■ h=0.01 
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